The calculation of multi-fractal properties of directed random walks on hierarchic 

trees with continuous branching 



(N 

O 

<N 

o 

fi: 

> 



O 

o 



m 
in 

(N 



X: 



David B. Saakian^'^'-^^ 

^Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan 
^ A.I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation, 
2 Alikhanian Brothers St., Yerevan 375036, Armenia and 
^National Center for Theoretical Sciences (North): Physics Division, 
National Taiwan University, Taipei 10617, Taiwan 

We consider tlie hierarchic tree Random Energy Model with continuous branching and calculate 
the moments of the corresponding partition function. We establish the multifractal properties of 
those moments. We derive formulas for the normal distribution of random variables, as well as for 
the general case. We compare our results for the moments of partition function with corresponding 
results of logarithmic 1-d REM and conjecture a specific powerlaw tail for the partition function 
distribution in the high-temperature phase. Our results establish a connection between reaction- 
diffusion equations and multi-scaling. 
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INTRODUCTION 



Random energy models on hierarchic trees are rather 
well-known and much investigated objects of statistical 
physics [H,[2|i|ll:[3l- The tree is constructed via a deter- 
ministic branching process, starting from the root node 
and adding q branches from every node of the tree, so 
that after K steps of the procedure there are nodes 
of the last generation (the endpoints of the tree). One 
then associates random energy variables to every branch 
of the tree, and further attributes the random energy to 
every node of the last generation by adding up all the en- 
ergy variables along the unique path connecting the given 
endpoint to the root. At the last step of the procedure 
one builds up the partition function [ij-Q. In the con- 
text of statistical physics the model has been identified 
with directed polymers on the trees jand has been also re- 
lated to the spin-glass model REM Independently, 
similar models have been applied and extensively inves- 
tigated also in the context of financial mathematics and 
turbulence , , HI] , [l^] , [H • They are also intimately 
connected to 2-d conformal models The hier- 

archical structure of the tree naturally induces recursive 
relations for the partition function, with the branching 
number q being a parameter in the recursion. While by 
construction the branching number q is an integer, one 
can formally consider the recursive equations in 
the g -> 1 limit, simultaneously allowing isT ^ 1 and 
keeping the number of endpoints q^ fixed. Such an idea 
has been suggested in [isj . later worked out in more de- 
tail in [l6| . Along that line an exact renormalization 
equation for continuous tree models has been derived in 
[l7| for the general distribution of random variables. 

The hierarchic tree models belong to the type of Ran- 



dom Energy Models (REM's) which all share the leading 
term in the free energy with that of the simplest Der- 
rida's model. The investigation of REM like models in 
finite dimensions has been started in [l8| and got serious 
impetus from recent solution of 1-d models with loga- 
rithmic correlation of energies at different sites [l^, [lo] . 
That solution essentially used the generalization of Sel- 
berg integrals [2l| , to complex number of integrations 
(see rigorous mathematical justification in [l^l). The Id 
models are intimately connected with the model [23|.[25j. 
Recently [2^ we have developed a statistical physics ap- 
proach to related dynamical MSM models ^3,- While 
all three types of REM models: hierarchical tree 1- 
d logarithmic REM Q and MSM [13 have exactly the 
same mean free energy as the standard REM, they have 
different distributions for the free energy and the parti- 
tion function. Most essentially, in the Derrida's REM the 
probability distribution of the partition function has no 
fat tails in the high temperature phase Q , while such tails 
are present in Id models [l9l . [20| . In this work we give 
indication that such tails exists in the hierarchic model, 
and also reveal the multi-fractal properties of the later, 
which is important for applications (7j-[To|. For a recent 
discussion of multifractal properties of REM, and other 
models with logarithmic correlations see [1^ [2^ and ref- 
erences therein. 

In multifractal approach (30j.[3l| one considers the mo- 
ments of a random variable (partition function) Z at some 
scale I defined with respect to the maximal scale L: 



(1) 
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where the exponent ^(n) defines the multi-scaling. 

Knowing the ^(n) and the coefficients Z„ 0{l) in 
the limit of large L, one can reconstruct the probability 
distribution of the partition function. In this way such a 
distribution has been explicitly derived for 1-d logarith- 
mic REMs [H, [13 

In the present article we will calculate ^{n) for hierar- 



2 



chical tree, and derive recursive relations to define 
which, in principle, could be applied to calculate the 
probability distribution of Z within some accuracy. We 
shall see indications of the fact that the divergence of 
moments is essentially the same as for lo? case, thus the 
two types of models must share the same fat tails. 



II. THE CALCULATION OF THE MOMENTS 

A. The model with Normal distribution of random 
variables 

Let us define the model outlined in the introduction 
following the papers Starting with the hierar- 

chical tree with integer branching number q and end- 
points, we consider two such end-points Wi and Wj. The 
two paths connecting these points to the root coincide 
up to a level m, counting from the root. Accordingly, we 
can define the hierarchic distance between the two points 




where Vq is defined as 



q'' ^ e'" 



mVp 
K 



L 



(2) 



(3) 



Associated with every branch of the tree is a Gaussian 
random variable e distributed according to the law 



K 



exp[ 
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(4) 



We define the energy yi at the endpoint Wi as a sum of 
corresponding variables e sampled along the unique path 
connecting the point Wi to the of the tree. 
We further define the partition function 



Z = Y^e 



where the sum is over 
Then we obtain 



(5) 



I end-points. 



< Z(e^,6)" > 



< e 



(6) 




FIG. 1: The sum is over three marked endpoints /, J, K of 
the hierarchic tree. Only the Fig. l.a contributes to the 
leading order. The locations of /, J, K given by Fig. l.b yield 
negligible 0{q — 1) contribution to the partition function in 
Eq.(8). 



We formally replace the sum over the tree with the 
integration over a measure dw [H, [l^ , 



< Z(e^,6)" 



n „ 
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> 
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We are able to calculate < Z" > for the integration range 
going over the maximal distance V. When calculating the 
correlations in Eq.(8) we use the following trick: the part 
of the trajectory of Wi with the hierarchic length v which 
has no overlap with trajectories from other points yields 
a factor while the part of the trajectory common 

to n such trajectories and hierarchic length v yields a 
factor e"''^'''/^^ 

Consider n ~ 2 for simplicity, with summation going 
up to the maximal distance V. We then have e^~^ po- 
sitions for the top level of hierarchy. Then, summing 
(or "integrating") over the number of two points at the 
distance e^" we get 

< Z{e^,bf >-- 



^„gV-t-'»g2(Vo-^)/3'+f/3=-Vo/3^ 



1 



^V{2+b)-bVo 



1 



(9) 



where we have introduced the parameter b = —0^. 

Consider now the integration in Eq. ([5]) going up to the 
maximal hierarchic distance V. We assume the following 
Ansatz: 

< Z(e^, 6)" >= e(^-^»)("+^"("-i)/2)e"^oZ„ (10) 



In Eq.([6|) the sum is over all endpoints of the tree sepa- 
rated by the hierarchic distance V from a chosen point. 
Now we consider the limit 



1 



(7) 



while keeping the value of Vq fixed and large. Recall that 
Vq = K\ogq. In such a limit there are e" end points at 
the hierarchic distance v. While calculating the < Z" >, 
we may discard the contributions given by Fig. lb. 



where we have introduced dimensionless coefficients Z„. 
The results of our calculation support the ansatz in Eq. 
(Uni), see Eqs.([13l)-(I14l). 

Actually the ansatz Eq. pH]) is correct only when 



(n-1) 



b 

n— > 
2 



(11) 



otherwise there is a phase transition to a different phase 
0, which manifests itself via diverging integration in 
Eq.®. 
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We will explicitly consider the case of positive b (imagi- 
nary /3) while calculating some eventual expression (there 
is no any restriction on n like Eq.(ll) ), like the analytical 
approximation of the fractional moment 



<Z''{V,b) >, 



(12) 



and then continue analytically the resulting expressions 
to the realistic case of negative 6 < 0. The analytical 
continuation gives wrong results in SG, when we con- 
tinue the expressions for moments to the other statistical 
physics phase. Hopefully in our case we arc interested to 
continue to expressions of moments in the same phase. 

Considering now the general case we calculate Zn re- 
cursively. At the highest hierarchy level, n points are split 
into two groups with m and (n — ni) points, 1 < m < n ac- 
cordingly, see Fig. 2. Having the expressions of Z^i^ Zj^^^ji 
at our disposal, we can calculate Z„. Wc assume the min- 
imum hierarchy level where all the paths from the root 
to n points meet is given by v. There are e^~" such lo- 
cations on the tree. Integrating over the positions of all 
points with a given v we arrive to 



l<m<n 



d^;[Z™(i>)Z„_™(i>)]e^-''+'"'"("-") (13) 

Jo 

Using the scaling Ansatz Eq. ([TOl) gives 

l<m<7i 



Zm{b)Zn-m 

71 — 1 + ^{m{m — 1) + (n — m){n — m — 1) — 2m{n — to)) 



El<m<«^'»(^)^"-m(&) 

{n - 1)(1 -^n5/2) 



(14) 



Thus, the problem amounts to solving the recursive equa- 
tion with the initial condition ((T^ . 

Eq. (fTO|) defines the so-called multifractal scaling. Wc 
can rewrite it in the form 

Let us now define a characteristic function 

n=l 

for those values of /3 where the sum converges. We rewrite 
Eq.(14) as 

(n-l)(l+n6/2)i„(&)= ^ Z,^{b)Z,,_^{b) (17) 

l<m<n 




FIG. 2: The graphics for calculation of Zn recursively frac- 
turing n = I + m. The groups with I and m endpoints meet 
at the hierarchic level v at the point A. 



Using the equations: 



l<m<n 



E7iZ„(6)x"=x— Ez„(6)x", 



E 



n'Z„(6)x" = (x— )' E 4(6)x", (18) 



we derive the following ODE for the characteristic func- 
tion: 



bx^ (Pu{x) du 2/ \ 



m(0) = l,it'(0) = 1 



(19) 



Unfortunately we could not solve Eq.(19) to derive ex- 
plicitly analytical expressions for Zn in case of general 
(complex) values of n, which precludes from following 
the procedure of [l^,[20|. Nevertheless, Eq.(19) allows us 
to extract a few first moments. We find: 



1,^2 



Z4 = 



1 ^ _ 2 
TT&' ^'^ " (l + &)(2 + 3&) 
1 



(l + 6)(2 + 3&)(3 + 65) (1 + 5)2(3 + 66) 



(20) 



The above relation allows us also to calculate the 
asymptotic of Zn at large n. In doing this we assume 
that our function u{x) has a singularity at some p of the 
form [H 



u{x) 



(1 _ £)a 



(21) 



Putting the latter expression into Eq. p^ . we derive: 



a = 2, 
362 

P 

Eg. ip^ gives the asymptotic expression 
^ _ 3(?i+l)62 



(22) 



(23) 



4 



We can define the value of p only numerically, looking at 
Zn for large values of the parameter n. A similar prob- 
lem has been considered in [sH . Slightly modifying their 
results we arrive at the following algorithm to calculate 
p. For a given integer ng we are looking the minimum 
over J and obtain: 



P 



6(n~l + n(n-l)bV2) 



0<i<no 



(24) 



It is possible to get accurate values of p by increasing the 
value of ng. 



In [15|,[17[ the following reaction-diffusion equation 
has been derived for the same model: 



dG{x,v) d^G{x,v) 



dv 



+ G{x,v)hiG{x,v) (25) 



where Q < v < V played the role of time parameter 
in the reaction-diffusion equation. Wc should solve that 
equation with the initial value of G 



G(x, 0) = exp[- 



(26) 



Therefore, there must be a relation between ODE p9|) 
and PDE (^5)) : solving Eg. ipsj) one can find the solu- 
tions of Ea. ([T5)) . Eq. ips)) allows one to be investigated 
via the traveling wave approach [l7| and fP =2\s its crit- 
ical point. The critical point of Eq.(19) should therefore 
correspond to h = —2, and is associated with the an- 
ticipated freezing transition from the high temperature 
phase to spin-glass-like phase. 



B. The case of general distribution 

Instead of Eq.(4) we now consider 

Vo 1 f°° Vq 

= ^ J dh cxp[—(j){ih) - ihx] (27) 

Then for the the sum yi of K random variables x, we can 
compose the distribution p(Vb,e) simply by multiplying 
(j) in the exponent of Eq. (27) by K: 

p(Vo,e) = ^ f dh eMVoHih) - the] (28) 



We define 



-iSvi-VoHm 



(29) 



We are calculating the correlations in < Z" > using 
the following trick: the part v of the trajectory of Wi, 
which belongs only him (no intersection with trajecto- 
ries of other points), gives a factor e'^^''^", while the part 
V of trajectory common to n trajectories, gives the factor 

0(ti/3)d 



Repeating the calculations of the previous section we 
obtain: 



<Z' >= 



^yg(V'-t-)</.(2/3)g2«0(/J)-2yo0(/3)g2t-gy-t;(^3Q) 



The factor eS^~^'>'^^^^^ corresponds to the part of trajec- 
tories [A, O], see Figl.a, while any of the lines [I, A] and 
[J, A] gives the factor e^^'^^^\ with e^" being the number 
of possible locations of /, J at the hierarchic distance w, 
whereas e^"" is the number of possible positions of the 
point A. 

Performing the integration over all < w < T^, we 
arrive at 

Zj(e'',« = z,(;3)e<''-''°>'=*<''>-*'"'+'> 

Let us derive recursive equations to calculate Zn,n > 1, 
defined as 

Z„(/3, V) = e(^-^o)M('3)-<#-M)+»)e"^oi„(/3) (32) 

We need to consider all possible splitting n ~ m+{n—m), 
with 1 < m < n. Let as assume that two group with m 
and n — m endpoints of our hierarchic tree are separated 
by the hierarchic distance u, and their trajectories meet 
at some point A. The calculations similar to those used 
to derive Eq. (fT4| give: 

Z„(/3,\/)= ZmiP 



1< ?n < n 



^(y_„)(l+0(„/3)_„0(/3)) 



Then the scaling in Eg. ([5^ yields: 

I]l<m<n Zm{P)Zn-m{l3) 



Zn{P) 



n - 1 - (j){nl3) + n(f){l3) 



(33) 



(34) 



One can use Eqs. (|5T|) . (|34l) to calculate any positive inte- 
ger moment of the partition function under the condition 

n0(^) + n < 1 + (t){np) (35) 

Eq. ([32|) implies the multifractal behavior with 

^n) - nm - 0(n/3) + 1 (36) 

For the corresponding generating function we now ob- 
tain an equation 



(l + 0(/3))x 



du{x) 
dx 



x-^)u 

dx 



u{x) + v?{x), 
u(0) = 1,m'(0) = 1 (37) 

For the same model in jlTj . another equation has been 
derived: 



dG{x,v) 
dv 



(l){dx)G{x,v) +G{x,v)\nG{x,v) (38) 



Thus two equations must be related. 
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C. Comparison of Logarithmic 1-d REM and 
hierarchic model with normal distribution 

In (20I I the following partition function has been con- 
sidered: 



Z = t^ I dye^ 
la 

\x — x'\ 



< y{x)y{x') >= 2 In ■ 



(39) 



with the parameter e being an ultraviolet cutoff needed 
to regularize the model. 

The above model shares the multifractal properties 
with our model for normal distribution of random en- 
ergies, with the mapping 



2/32 = /32 



(40) 



Let us compare the moments. Ref.[20| gives the following 
expression for the moments: 



n 



r(i - Q- - i)7)"r(i - 7j) 

L r(2-(n + j-2)7)r(l-7) 



(41) 



where 7 = /S^. 

We have, naturally Zi = 1 and further 



r(i - 27) 



1 



r(2-7)r(2-27) (l-7)(l-27) 

2 r(i~27)3r(i-37) 

^ r(2 - 27)r(2 - 37)r(2 - 27) 

r(i - 27)3r(i - 37)3r(i - 47) 
r(i - 7)r(2 - 37)r(2 - 47)r(2 - 57) 



(42) 



Comparing with the results of the section II-A, we sec 
that second moment Z2 has different expansion for the 
small j3'^, 



Z. 



1 + 37, 



(43) 



More important information contained in moments is 
however the smallest real pole 7„ as a function of 7 (re- 
spectively, h). Indeed, precisely those poles define the 

— 1/2 

critical temperatures r„ 7„ below which the given 
moment of the partition function start to diverge. It 
is easy to see from Eq. (|4T|) that for one-dimensional 
model 7 — 1/n, n ~ 2,3... while Ea. (|l^ gives the 
smallest pole poles at 6/2 = 1/2 for Z2, b/2 = 1/3 
for Za, h/2 = 1/4 for Z4, etc. The results for a few 
lower moments indicate that the two models share the 
same sequence of " transition temperatures" . Though 
we are unable to prove this statement in full general- 
ity, the factor (1 — bn/2) appearing in the denominator 
of the right-hand side of the recursion Eq.([T4]) makes the 
statement very plausible, if no special cancelations hap- 
pen. Assuming that correspondence is correct, the lat- 
ter property would imply the probability density for the 



partition function in both models share the same power 
—1— ^ 

law tail: ViZ) (x Z ^ valid everywhere in the high- 
temperature phase /3 < /3c = 1- Such a tail was indeed 
proposed as one of the universal features characterizing 
the class of models with logarithmic correlations [20j . 

We obtained analytical expressions for the hierarchic 
model's Z„ at positive b, Eq.(14). Having large series of 
Zn for positive b, one can try to construct an approximate 
expression for the probability distribution and fractional 
moments. 



III. CONCLUSION 

The investigation of the statistical mechanics of REM- 
like models in finite dimension and their multifractal 
properties remains an active field of research, especially 
after the advances achieved in 27 1 and In this ar- 
ticle we confirmed that the multifractal properties are 
shared by the third class of REM-like models: directed 
polymers on a hierarchic tree. We calculated the cor- 
responding moments important for the applications. In 
particular, we use those moments to conjecture the power 
law tail of the distribution of partition function Z in the 
high-temperature phase. Unfortunately it is impossible 
to derive exact probability distribution, as has been done 
in for the 1-d logarithmic REM. Nevertheless, our 
formulas allow the calculation of infinite series of the mo- 
ments for some values of parameters. We hope that such 
information could be applied to recover the probability 
distribution, a well known problem in probability theory 

Our results are interesting for the mathematics of 
reaction-diffusion equations as well, as we related them to 
a certain nonlinear ODE. Moreover, we found the connec- 
tion between the multifractal spectrum and the hierar- 
chic tree model Eq. (pS)) , while the latter has been mapped 
to some reaction-diffusion equations [l3, see Eq. ([55)) in 
the present article. Thus the reaction -diffusion equa- 
tions have some multiscaling structures. One can con- 
nect reaction diffusion equation and nonlinear ODE with 
the interesting versions of multi-fractal phenomenon, i.e. 
[37| . and try to investigate them for deriving the multi- 
fractal scaling from the reaction diffusion equation (|38|) . 
Such a work is currently in the progress. 

The directed walk models on hierarchic tree describes 
a rather rich physics, as well as some connections with 
quantum models in finite dimension. The relation to 2-d 
conformal field models is rather well known [l3i[l3i and 
the work [l3| mentions the relations to quantum disorder 
problem in finite dimensions [36j . It will be interesting to 
try to apply the methods of the current work or our renor- 
malization group psp to the finite dimensional quantum 
disorder problem. 

I thank Y. Fyodorov and V. Zacharovas for the discus- 
sion, Acadcmia Sinica for support. 
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